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Abstract 

We compute the Schrodinger functional (SF) for the case of 
pure SU(3) gauge theory at two-loop order in lattice perturba- 
tion theory. This allows us to extract the three-loop /3-function 
in the SF-scheme. These results are required to compute the run- 
ning coupling, the A-parameter and quark masses by finite size 
techniques with negligible systematic errors. In addition, we can 
now implement two-loop 0(a) improvement in simulations and 
extend and study series in alternative ( "tadpole-improved" ) bare 
couplings. 



1 Introduction 



Traditionally, in computer simulations of lattice QCD physicists have been 
concerned with the computation of low energy quantities which are inacces- 
sible to perturbation theory. Typical examples are ratios of hadron masses, 
the static quark potential at large separations and glueball masses. At high 
energy, where QCD was first established as a promising theory of strong in- 
teractions, perturbation theory in the dimensionally regularized continuum 
theory is the standard computational tool. Here it provides predictive power 
by reproducing a host of high energy data in terms of a renormalized run- 
ning coupling and running quark masses at that energy scale. As a result the 
coupling in the MS-scheme at the scale of the mass of the Z-boson is fixed to 
be oty^(M z ) ~ 0.11. For a long time the two energy sectors — though they 
are believed to belong to one and the same theory — have coexisted rather 
independently. Only in recent years was the task of connecting them by 
computing in the theory that is completely parameterized at low energy, 
taken up (see |l|, ||, || for reviews). 

Referring also to low energy this result has to be obtained by simulation 
on the lattice. The extraordinary difficulty of this computation is caused 
by the fact that we seem to have to simultaneously control the continuum 
limits of two widely separated physical energies. Naively this requires un- 
feasibly large lattices. The problem is circumvented by the somewhat indi- 
rect strategy developed by the ALPHA-collaboration in recent years @, |3|. 
Here we numerically trace the evolution of the coupling «sf defined by the 
Schrodinger functional from low to high energy. While not being a directly 
physical quantity, it is nevertheless expected to possess a well defined con- 
tinuum limit and the motivation for this particular choice is explained in 
||. Due to its definition on a finite system, it runs with the size L. Once 
its evolution is known, it remains to connect both its low and high energy 
values to more experimentally accessible quantities. At low energy (i.e. at 
large physical volume) this is done by numerical simulation, which however 
no longer involves unmanageable scale ratios. At high energy, asF on small 
systems can be perturbatively related to a^fg- It has turned out that with 
the precision that is both experimentally required and achievable in the var- 
ious simulation steps, this connection is needed at two-loop order. Due to 
the difficult structure of lattice perturbation theory one usually has to con- 
tent oneself with one- loop results. In fact, the results presented here and its 
analogue for gauge group SU(2) are the only lattice two-loop results of 
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renormalized quantities known to us. 

The perturbative relation between cksf and is established in two steps 
through the bare lattice coupling go as an intermediate quantity. The two- 
loop relation between and go has been worked out in || for Yang-Mills 
theory with gauge group SU(iV) for general N. The present paper completes 
the two loop relation for quenched QCD by deriving the connection between 
g Q and asF for the realistic case of SU(3)[|. The analogous SU(2) calcula- 
tion was reported in ||. We note that here Feynman diagrams have to be 
evaluated and extrapolated numerically. This cannot be done for general N, 
and the more complex group structure of SU(3) represents a very significant 
complication. 

After this introduction, we define the problem for SU(3) and give a few 
hints to details of the computation in section 2. In section 3 we compile 
results for finite lattices and analyze and extrapolate them. Finally in section 
4 some applications of our results are reviewed. Some preliminary results of 
this work were already given in ||. 



2 Perturbation theory for gauge group SU(3) 

we compute the free energy T 



In the Schrodinger functional approach 
given by the lattice path integral 



D[U]e 



-5(10 



(2.1) 



The geometry is a finite box of size L with periodic boundary conditions 
in the spatial directions and Dirichlet boundary conditions in time. In the 
lattice realization the latter means that the spatial links U(x, k),k = 1, 2, 3, 
are frozen to fixed SU(3)-elements for the layers at time coordinate x° = 0, L, 



U(x, k)\ x o =0 
U(x,k)\ x0=L 



la 



exp(aC) = exp — diag( 







exp(aC') = exp — diag 



HI 



'1j V2> V3, 



(2.2) 
(2.3) 



Here a is the lattice spacing, and 4>u 4>'i are certain [|J dimensionless numbers 
that depend on one free parameter 7] that we shall vary later on to probe the 
system. Both our perturbative calculation here and simulations reported in 



1 Results reported here are based on ref. where more details can be found. 
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and more recently in (TO, ITTJ are restricted to the choice "A" of [ffl for 



these numbers. All other links interior to the box are integrated over with 
the invariant SU(3) measure in ( |2.1[ ). The action is defined by the usual sum 
over oriented plaquettes, 

S(U) = ±^w(p)tr(l-U p ). (2.4) 
9o P 

The weight w{p) is unity for all plaquettes except those at the boundary that 
contain the time-direction and one of the frozen spatial links where we put 

w( P ) = c t (g ) = 1 + # ] gl + cf ] g 4 Q + ... (2.5) 

The freedom of adjusting this weight is required for improvement of 0(a) 
lattice artefacts that are otherwise introduced by the surfaces. 

We set asF = ^? 2 /4vr and the SF-coupling g is defined from the response 
in the energy T to infinitesimal changes in the surface fields by varying rj, 

f = k/V , (2.6) 

where T' is the derivative with respect to i] at rj = and A; is a constant. It 
is fixed by normalizing the leading term in the perturbative expansion 

f{L) = g 2 Q + m^L/^g 4 + m 2 (L/a)g 6 + ... (2.7) 

The one loop coefficient mi was computed in ref. |4[]. Our objective here 
is the computation of m,2 for gauge group SU(3). In order to determine the 
improvement coefficients in ( |2.5| ) we need to explicitly exhibit the dependence 
of the coefficients m; on them. Inspection of the structure of the contributions 
leads us to write 



mi 



m\ + c t (1) m5, (2.8) 



m 2 - mf = + c\ L, m b 2 + c f (1) " m\ + ci ] m d 2 . (2.9) 

To work out the perturbative expansion of T we fix the gauge as discussed 
in ||. It is an important advantage of the SF-framework that there is a 
unique background field of minimal action that interpolates between the 
surface values 



V(x, k) = exp[a(C + (C - C)x°/L)), V(x, 0) = 1, (2.10) 
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and all fluctuation modes receive gaussian damping. To deal with the SU(3) 
Lie-algebraic structure of the fluctuations it is advantageous to introduce a 
suitable basis. We used generators proportional to 

T° = diag(0,l,-l) oc (C'-C), 
T 1 = diag(2,-l,-l), (2.11) 
(T<*% = <^f; (a(3) = (12), (13), (23), 

such that T°, T 1 span the Cartan subalgebra. The T a @ and their adjoints 
are natural generalizations of the ladder-operators in SU(2). Now covariant 
derivatives in the background field acting in the adjoint representation, — an 
ubiquitous structure in our perturbation theory — are diagonalized in group 

space, 

[T°,T al3 } = h af3 T a/3 (2.12) 

with numbers h a p. In this basis the numerical construction of the propagators 
in the background field is reduced to the SU(2) case. We just have to apply 
the procedure described in section 5 of ref. [||] several times for various 
effective background field parameters. 

While the vertices for SU(2) were still 'hand-coded' in Fortran, the SU(3) 
structure was so much more complicated that we let Maple do the work. 
The three-gluon vertex, for instance, consists of 532 terms. As for SU(2), we 
performed the summations in position space, so that the most demanding big- 
mac diagram shown in fig. [1] has only 0((L/a) 5 ) instead of 0((L/a) 8 ) terms 
in momentum space. Its total number of contributions for L/a = 32 after 




Figure 1: The most computer-time intensive diagram, the big-mac with 
three gluon propagators. 

all possible reductions is however still about 5 • 10 12 . Along with each vertex 
amplitude and propagator we computed their exact ^-derivatives to form T'. 
To enhance the confidence in our numbers we took the same precautions as 
in ||: test of all symmetries before using them to reduce the number of terms 
by a factor of 36 altogether, check of gauge independence and comparison 
with a slow but independent code for the smallest lattices. 
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3 Extraction of results 



Some of the coefficients in ( |2.9| ) can be given in closed form. Since the 
background field ( |2.10|) leads to constant plaquettes it is easy to see that 



m\ = mi = ~. (3.1) 

The tree-level contribution m 2 can be evaluated in closed form. The en- 
suing expression is however a complicated function of the background field 
parameters and for our purposes the leading behavior is sufficient 

m c 2 = |-^ + 0(a 5 ). (3.2) 

The remaining contributions m", m 2 and m 2 are tabulated in table |l[ The 
generically expected behavior for 2-loop lattice Feynman diagrams suggests 
an asymptotic expansion 

oo 

™l = + SnHL/a) + t n \n 2 (L/a)}(a/L) n . (3.3) 

n=0 

For m} we expect the same structure without ln 2 -terms and the same is true 
for the two-loop diagrams contributing to m 2 which in addition starts at 
n = 1 only. The coefficients of these expansions are extracted from the series 
of finite lattices and errors are estimated by the blocking procedure described 
in ||12| . The Callan-Symanzik equation 

L -^?{L) = ~m = bof + hf + b 2 g 7 + . . . (3.4) 

implies to = and fixes all logarithmic divergences in terms of the universal 
coefficients 

11 102 
b ° ~ (4^' 6l " (4^- (3 ' 5) 

These values were first verified to high accuracy from the data and then 
assumed to be exact to extract subleading terms. Cancellation of all terms 
of 0(a) in mi and m 2 is achieved to our numerical accuracy by choosing 

Ct = 1 - 0.08896(23)^ - 0.0301(25)^ + O(^). (3.6) 
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L/a 


m\ 


m\ 


i, 

m 2 


4 





529761318876020 


0.052943274533542 


0.106270501583520 


5 





564680831258396 


0.056838593961882 


0.100438272480780 


6 





593212325999814 


0.05975069070789 


0.09439258604949 


7 





617378117870096 


0.06197806027570 


0.08867849757097 


8 





638244418293146 


0.06374136571959 


0.08343000201173 


9 





656541534280898 


0.06518532139185 


0.07867547537704 


10 





672800386347018 


0.06640139951338 


0.07439325408456 


11 





687413016523170 


0.06744855730398 


0.07054052556373 


12 





700673174318993 


0.06836628672255 


0.06706893587917 


13 





712804791119344 


0.0691821005300 


0.06393184996835 


14 





723981656094021 


0.0699158333750 


0.06108717717361 


15 





734340962964557 


0.0705822007834 


0.05849809311936 


16 





743992719364831 


0.0711923868147 


0.0561328521011 


17 





753026381769789 


0.0717550680760 


0.0539642582909 


18 





761515615421704 


0.0722770966172 


0.0519690528035 


19 





769521768768977 


0.0727639684486 


0.0501273255443 


20 





777096451745942 


0.0732201529077 


0.0484219921318 


21 





784283478608195 


0.0736493291767 


0.0468383453838 


22 





791120352757074 


0.0740545593455 


0.0453636778366 


23 





797639416324361 


0.0744384171842 


0.0439869669949 


24 





803868750843028 


n r\T a onono c a i on 

0.0748030854189 


0.0426986139530 


25 





809832890636312 


0.0751504302247 


0.0414902265431 


26 





815553393561204 


0.075482058988 


0.0403544392531 


27 





821049301867167 


0.075799365607 


0.0392847633581 


28 





826337517515302 


0.076103566396 


0.0382754618367 


29 





831433110259494 


0.076395728812 


0.0373214446266 


30 





836349572396415 


0.076676794651 


0.0364181805963 


31 





841099030858248 


0.076947598928 


0.0355616232918 


32 





845692424917554 


0.077208885368 


0.0347481480648 



Table 1: List of L/a-dependent coefficients m", and m b 2 . Numerical errors 
are beyond the quoted digits. 

With these values, referred to as 2-loop perturbative O(a) improvement, we 
finally quote 

mi = 26 ln(L/a) + 0.36828215(13) + 0(a 2 ), (3.7) 
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m 2 -m\ = 26iln(L/a) + 0.048085(63) + 0(a 2 ), (3.8) 

for the relation (|2.7| ) between g\ and g 2 . 

The two-loop connection between two different couplings determines the 
difference between the non-universal three-loop coefficients of their respective 
/5-functions. Hence we can employ the results in |§, where use is made of 
the three-loop /3-function in the MS-scheme to determine 

0.4828(88) 

62 = (3 ' 9) 

for the third coefficient in (|3.4| ). This is of the order of magnitude of the 
effective two-loop coefficient 63 = l-5(8)/(47r) 3 which was used in ref. [ffl to 
represent the data. 

A central quantity in the ALPHA-collaboration's approach is the step 
scaling function which generalizes the /3-function to finite rescalings, 

a(s,u)=g 2 (sL)\ u= - 2{L) . (3.10) 

On the lattice, it emerges as the continuum limit of a finite lattice spacing 
approximant E, 

a(s, u) = lim E(s, u, a/L). (3-11) 

All perturbative information about the convergence speed for s = 2 is con- 
veniently summarized in coefficients 81 and 8% defined by 

8^ (u, a/L) = n^a/L)-a{2,u) = ^ ^ + ^ ^ + _ 

(t(2. u) 

(3.12) 

where the upper index k — 0, 1, 2 refers to varying the degree of improvement 
by using an action with one, two or three terms in the expansion of q in 
Our results are collected in table 0. 



4 Applications 

The perturbative coefficients, whose computation we described in the pre- 
vious sections, have already found some applications which we now briefly 
discuss. Once the SF-coupling g is known for small box-size, it is to be con- 
verted to the MS-coupling at high energy. One could conventionally choose 



7 



L/a 






8? 




r(0) 
°2 






v(2) 

°2 


4 





01192 


-0 


01023 





00642 





00577 


-0 


00168 


5 





01155 


-0 


00617 





00514 





00501 


-0 


00095 


6 





01089 


-0 


00387 





00412 





00435 


-0 


00061 


7 





01004 


-0 


00262 





00334 





00382 


-0 


00043 


8 





00918 


-0 


00189 





00275 





00341 


-0 


00032 


9 





00841 


-0 


00144 





00230 





00307 


-0 


00024 


10 





00773 


-0 


00113 





00195 





00279 


-0 


00019 


11 





00714 


-0 


00091 





00168 





00256 


-0 


00015 


12 





00663 


-0 


00075 





00145 





00236 


-0 


00012 


13 





00618 


-0 


00063 





00126 





00219 


-0 


00010 


14 





00579 


-0 


00054 





00111 





00205 


-0 


00008 


15 





00544 


-0 


00046 





00098 





00192 


-0 


00007 


16 





00513 


-0 


00040 





00087 





00180 


-0 


00006 



Table 2: The perturbative lattice effects. 



the mass of the neutral weak boson Mz as a scale here. Another choice - 
attractive for asymptotically free theories — is to extract the A-parameter 
which is simply related to the behavior at asymptotically large energy. It is 
a renormalization group invariant given by 



A SF = L-\b fr b ^e- 1 ^ b0 ^ x exp j- J° dg 



1 1 h 



(4.1) 

The conversion to A^ then amounts to an additional known factor. If we 
insert a small g(L) belonging to a very small L into formula (|4.1| ), then the 
exponentiated integral is close to unity. In the A-parameter in the MS 



scheme at zero flavor number has been estimated to be 

A|| = 251(21)MeV. (4.2) 

The quoted error is the statistical error from the simulations of the couplings. 
The three-loop term of the /3-function has been used and this renders the 
systematic error due to further orders negligible. Without the three loop 
term, a hard to quantify systematic error of the order of the statistical one 
would have remained. 
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In (T3| the continuum limit of L maiX /r = 0.718(16) is computed. The box- 



size L max is defined [|J via the Schrodinger functional by g 2 (Z/ max ) = 3.48. 
The distance ro is fixed though the relation t^F^q) = 1.65, where F is the 
force between static quarks |14]| . It is expected to be close 0.5 fm in physical 



units. Some of the lattice artefacts derive from the definition of L max and are 
expected to be 0(a). Series of simulations have been conducted both with 
the one-loop and two-loop forms of c t , eq. (p.6|). In fig. 5 of |13| the effect 
of Cf 2 ^ is very clearly visible. It practically eliminates the 0(a) term in the 
continuum extrapolation which thus becomes more reliable and accurate. 

We conclude with some remarks on the numerical quality of perturbative 
expansions in the bare coupling, ao = go/^, or improved bare couplings 



| 15fl like dp = a /P, where P is the average plaquette (in infinite volume) 
for the corresponding bare coupling. There are few occasions to investigate 
such series beyond one- loop order as we can do now in the case of agp. The 
well-known problem with bare perturbation theory shows up in the size of 
the coefficients at scale q = 1/a, 

a SF (l/a) = a + 4.628 a 2 + 29.0 ajj + 0(a£). (4.3) 

If we choose a scale s/a where the 1-loop term is cancelled, then s is unnat- 
urally large and the next term is still not small 

a SF (14.06/a) = a + 4.18 a 3 + 0(a£). (4.4) 

If we change to the boosted coupling the analogous formulas look somewhat 
better, 

a gF (l/o) = a P + 0.439 c4 + 2.43 c4 + 0(c4), (4.5) 
a SF (1.285/a) = a P + 1.91 a% + 0(a P ). (4.6) 

The same tendency is obtained if we minimise the two-loop coefficient instead 
of cancelling the one-loop term because the scale ratio r between the former 
and latter method comes out to be 

r = exp(6 1 /(462)) (4.7) 

which is close to unity. With typical couplings like «p ~ 0.11 for f3 = 6.5 it 
still is difficult to be confident about systematic errors incurred by truncating 
the series. 
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As the ALPHA collaboration moves to the inclusion of two massless fla- 
vors of quarks into simulations of «sf we plan to extend the present two-loop 
calculation by the required diagrams with quarks. We hope to report on this 
in the near future. 
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